Production of the X(3872) in B Meson Decay 
by the Coalescence of Charm Mesons 



Eric Braaten and Masaoki Kusunoki 
Physics Department, Ohio State University, Columbus, Ohio 43210, USA 



o 
o 

(N 

as: 
< 

m 

(N 

(N 
> 

o 
o 

a: 
i 

^ ■ 



X 



Shmuel Nussinov 

School of Physics and Astrononmy, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel 

(Dated: February 1, 2008) 

If the recently-discovered charmonium state X(3872) is a loosely-bound S-wave molecule of the 
charm mesons D°D*° or D*°D°, it can be produced in B meson decay by the coalescence of charm 
mesons. If this coalescence mechanism dominates, the ratio of the differential rate for B + — > 



D°D*°K+ near the D°D*° threshold and the rate for B 4 



XK + is a function of the D°D*° 



invariant mass and hadron masses only. The identification of the X(3872) as a D° D*° / D*° D° 
molecule can be confirmed by observing an enhancement in the D°D*° invariant mass distribution 
near the threshold. An estimate of the branching fraction for B + — > XK + is consistent with 
observations if X has quantum numbers J PC — 1 ++ and if J/ip tv + ty~ is one of its major decay 
modes. 



PACS numbers: 12.38.-t, 12.38.Bx, 13.20.Gd, 14.40.Gx 

The recent unexpected discovery of a narrow char- 
monium resonance near 3.87 GeV challenges our un- 
derstanding of heavy quarks and QCD. This mysterious 
state A"(3872) was discovered by the Belle collaboration 
in electron-positron collisions through the B-meson de- 
cay S ± — >■ K ± X followed by the decay X — > J/ip tt + it~ 
0. The discovery was confirmed by the CDF collabo- 
ration using proton-antiproton collisions 0]. The X is 
much narrower than all other charmonium states above 
the threshold for decay into a pair of charm mesons. Its 
mass is also extremely close to the threshold for decay 
into the charmed mesons D°D*° or D*°D°. 

The proposed interpretations of the X(3872) include 
a D-wave charmonium state with quantum numbers 
J = 2 or 2 h , an excited P-wave charmonium 
state with J = 1 ++ or l" 1 , a "hybrid charmonium" 
state in which a gluonic mode has been excited, and a 
D a D*°/D*°D molecule 0, US H0,& HM 111 El 
The possibility that charm mesons mi ght form molecular 
states was considered some time ago [lj, [l5j, [l6j . If the 
binding is due to pion exchange, the most favorable chan- 
nels are S-wave with quantum numbers J PC = 1 ++ or 
P-wave with h Q. The proximity of the mass of X to 
the D°D*° threshold indicates that it is extremely loosely 
bound. If X is an S-wave D D*°/D*°D molecule, the 
tiny binding energy introduces a new length scale, the 
D°D*° scattering length a, that is much larger than other 
QCD length scales. As a consequence, certain properties 
of the X I D° D*° I D*° D° system are determined by a and 
are insensitive to the shorter distance scales of QCD. This 
phenomenon is called low- energy universality. 

A challenge for any interpretation of the AT (3872) is to 
explain its production rate. This could be problematic 



for the identification of X as an S-wave D D*°/D°D* 
molecule, because it can readily dissociate due to its 
tiny binding energy. One way to produce X is to pro- 
duce D° and D*° with small enough relative momentum 
that they can coalesce into X. An example is the decay 
T(4S f ) — * Xhh', where h and h! are light hadrons, which 
can proceed through the coalescense into X of charm 
mesons from the 2-body decays of a virtual B and a vir- 
tual B. Remarkably, low-energy universality determines 
the decay rate for this process in terms of hadron masses 
and the width of the B meson ^3 ■ Unfortunately, the 
rate is supressed by a factor of (r b /ma) 2 and is many 
orders of magnitude too small to be observed. 

In this paper, we apply low-energy universality to 
the discovery mode B + — * XK + and to the process 
B+ -> D°D*°K+. We point out that the interpreta- 
tion of X as an S-wave D D*°/D*°D molecule can be 
confirmed by observing a peak in the D°D*° invariant 
mass distribution near the D°D*° threshold in the decay 
B + — > D a D*° K + . We also estimate the branching frac- 
tion for B + — » XK + . The estimate is compatible with 
observations if X has quantum numbers J PC = 1 ++ and 
if J/rp tt + 7t~ is one of its major decay modes. 

The mass of the X has been measured to be mx = 
3872.0 ± 0.6 ± 0.5 MeV by Belle 0] and 3871.4±0.7±0.4 
MeV by CDF 2]. It is extremely close to the D°D*° 
threshold 3871.2 ± 0.7 MeV. The binding energy is E b = 
—0.5 ± 0.9 MeV. If the state is bound, Eb is positive, so 
it is likely to be less than 0.4 MeV. This is the smallest 
binding energy of any S-wave two-hadron bound state. 
The next smallest is the deuteron, a proton-neutron state 
with binding energy 2.2 MeV. For two hadrons whose 
low-energy interactions are mediated by pion exchange, 



2 




FIG 

way. 



1: Feynman diagram for B + — * XK + via the first path- 



the natural scale for the binding energy of a molecule is 
m 2 /(2/i), where /z is the reduced mass of the two hadrons. 
For a D°D*° molecule, this scale is about 10 MeV, so E b 
is at least an order of magnitude smaller than the natural 
low-energy scale. 

If the binding energy of X is so small, low-energy uni- 
versality implies that the X / D° D*° / D*° D° system has 
properties that are determined by the D°D*° scattering 
length a and are insensitive to the shorter distance scales 
of QCD. The universal binding energy of the molecule is 



Eb = m D + m D * — m x ^ (2{ia 2 ) 1 , 



(1) 



where fi = mump* / (mo + ran*) is the reduced mass of 
the D° and D*°. The universal normalized momentum- 
space wavefunction at relative momentum k <C m^, 



iP(k) ~ (87r/a) 1/2 (fc 2 + l/a 2 Y 



(2) 



was used by Voloshin to calculate the momentum distri- 
butions for the decays X -> D°D°ir a and X -> D°D°j 
6]. The universal D°D*° elastic scattering amplitude at 
relative momentum fc cm <C is 



A[D°D*° -► D°D* } 



87rmi)m£). 
(i(-l/a- ik cm )' 



(3) 



where k CTn w [2/i(-E — — m^.)] 1 / 2 and £ is the total 
energy in the center-of-momentum frame. The amplitude 
A[D*°D° -> Z)°D*°] for scattering to the CP conjugate 
state differs by the charge conjugation C = ± of the 
channel with the large scattering length. Another conse- 
quence of low-energy universality is that as the binding 
energy Eb decreases, the probabilities for components of 
the wavefunction other than D°D*° and D* Q D Q decrease 



as E, 



1/2 



In the limit Eb — ► 0, the state becomes 



{\D*°D } ± \D°D*°))/V2 if C = ±. The rates for decays 
that do not correspond to the decay of a constituent D*° 
or D*° also decrease as E^ 2 . This suppression may ex- 
plain the surprisingly narrow width of the X. 

The decay B + — > AAT + proceeds through the weak de- 
cay b — > ccs at very short distances. The subsequent for- 
mation of XK + is a QCD process that involves momenta 



k as low as 1/a. The contributions from k ~ 1/a are con- 
strained by low-energy universality, but those from k > 
TOtt involve the full complications of low-energy QCD. 
We analyze the decay B + — » AAT+ by separating short- 
distance effects involving k > m T from long-distance ef- 
fects involving k ~ 1/a. The decay can proceed via the 
short-distance 3-body decay £?+ — > D°D*°K + followed 
by the long-distance coalescence process D°D*° — > X. 
It can also proceed through a second pathway consist- 
ing of 5+ -► D*°D°K + followed by D°D*° -► X. The 
amplitude for the first pathway can be expressed as 
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D°D* yl K 



0; 



xD(q + £,m D )D(q* - £,m D .) A[D°D*° X], (4) 



where g = {niD/'mx)Q and q* = (mp* /mx)Q are 4- 
momenta that add up to the 4-momentum Q of A and 
D(p,m) = (p 2 — m 2 + ze) -1 . The sum is over the spin 
states of the D*°. This amplitude can be represented by 
the Feynman diagram with meson lines shown in Fig. ^ 
We constrain the loop integral to the small-momentum 
region by imposing a cutoff \£\ < A in the rest frame of 
the virtual D° and D*°. The natural scale for the cutoff 
is A ~ m„. The amplitude for D°D*° to coalesce into X 
is determined by the D°D*° scattering length a: 



A[D°D*° -> X] 



= (l6TrZmxrnjjmE>* / /J 2 ^ 



1/2 



(5) 



where ex and e are the polarization vectors of X and 
D*° and Z is the probability for the X to be in a 
D a D*°/D*°D state. At the D°D*° threshold, the am- 



plitude for B+ -> D°D*°K^ 
invariance to have the form 



is constrained by Lorentz 



A[B+ -> D U D* U K+] = ci P ■ e 



(6) 



where P is the 4-momentum of the B meson and c\ is 
a constant. The amplitude for B+ -> D*°D°K+ has 
the same form with ci replaced by a constant Ca- In the 
D a D*° rest frame, the integral over £ of the two prop- 
agators in is proportional to the momentum-space 
wavefunction of X. The subsequent integral over £ is 
linear in the ultraviolet cutoff A for the low-momentum 



region: 



d i £ 



D(q + £,m D )D{q' - £,m D .) = 



ijiK 



•(7) 



(27r) 4 " 47r 2 mDmD» 

The total amplitude from the two pathways is 

A[B + -> AA+] = - (Zm x /7r 3 mz5?7i D -a) 1/2 

x( Cl ±c 2 )AP-e x . (8) 

The sign ± corresponds to the charge conjugation C = ± 
of A . Heavy-quark spin symmetry implies c\ — ci up to 



FIG. 2: Feynman diagram for B 
pathway. 



DD*°K + via the first 



corrections suppresssed by a factor Aqcd The in- 
terference is constructive if C = + and destructive if 
C = — . The dependence of the loop amplitude iJSJl on A 
is cancelled by a tree diagram with a B — XK contact 
interaction whose coefficient therefore depends linearly 
on A. If the X is predominantly a DD* molecule, there 
must be some value Ai of the ultraviolet cutoff for which 
the loop amplitude dominates over the tree amplitude. 
Squaring the amplitude, summing over spins, and inte- 
grating over phase space, the final result for the decay 
rate is 



XK' 



ZX 3/2 (m B ,m x ,m K ) 



64tt 4 



m B m x \ia 



Cl ±c 2 | 2 A2,(9) 



where X(x, y, z) = x 4 + y 4 + z 4 — 2(x 2 y 2 + y 2 z 2 + z 2 x 2 ). 

1 li 

Due to the factor 1/a, the decay rate scales like E h as 
E b -» 0. 

If another hadronic state H is close enough to the 
D°D*° threshold that X has a nonnegligible probability 
Zh of being in the state H, the decay can also proceed 
through a short-distance 2-body decay B + — > HK + . In 



1/2 
H 



this case, there is an additional term A[B + — > HK + ]Z 
in (JSJ. Its contribution to the decay rate also scales like 
e\^ 2 as Eb — * 0, because Zh scales like E^ 2 If 
C = +, one possibility for such a state is the excited 
P-wave charmonium state x c i(2P). Recent coupled- 
channel calculations of the charmonium spectrum suggest 
that (2P) is likely to be well above the D°D*° thresh- 
old We will henceforth assume that D a D*°/D*°D 
is the only important component of the wavefunction and 
set Z w 1. 

We can calculate the differential decay rate for B + — > 
D°D*°K + in the same way. There are again two path- 
ways: the short-distance decay B+ -> D°D*°K+ fol- 
lowed by the long-distance scattering D°D*° — > D°D*° 
and -> D*°D°K+ followed by D*°D° D a D*°. 
The amplitude for the first pathway can be represented 
by the Feynman diagram with meson lines shown in 
Fig. [3 The calculation of the amplitude is similar to 
that for B + — > XK + except that it involves the scatter- 
ing amplitude © instead of the coalescence amplitude 



■ m D - m D . (MeV) 



FIG. 3: The D°D*° invariant mass distribution for B + -» 
D°D*°K + for three different values of the binding energy of 
X. The distributions are normalized to 1 at fc cm = m^. 



@. In the loop amplitude for B + D D*°K + , we keep 
only the term © that is nonzero at the D°D*° threshold. 
There must be some value A 2 of the ultraviolet cutoff for 
which the loop amplitude dominates over the tree ampli- 
tude. The factor ci ± c 2 cancels in the ratio between the 
amplitudes for B+ -> D°D*°K + and B+ -► XK+. Our 
final expression for the differential decay rate is 



clT 



dM 



\B^ 



D°D*°K A 



DD' 



T[B + -> XK + ] 



A 2 



jJLO?k c 



A{ n(l + a 2 k 2 m ) 



where Mq d , is the D°D*° invariant mass and k c 
relative momentum in the D°D*° rest frame: 



(10) 



is the 



fc cm - \ 1 ' 2 {M DD , , m D ,m D ,)/(2M DD ,). (11) 

In (|1U[) . we have neglected terms suppressed by fc^ m /m|j. 
The invariant mass distribution is illustrated in Fig.|21for 
several values of the binding energy Ef,. The distributions 
are normalized to 1 at fc cm = m w . As the binding energy 

— 1/2 

is tuned toward 0, the peak value scales like E. ' and 
the position of the peak in Mq d » — (mo + mo*) scales 
like Eb- The observation of such an enhancement near 
the D°D*° threshold would confirm the interpretation of 
AT as a D°D* a /D* a D° molecule. 

The Babar collaboration has recently measured the 
branching fractions for B + to decay into D°D°K + , 
D°D*°K + , D*°D°K+ 1 and D*°D*°K+ to be (0.19 ± 
0.03)%, (0.47±0.07)%, (0.18±0.07)%, and (0.53±0.11)%, 
respectively 0. We use these measurements to esti- 
mate the branching fraction for B + — > XK + . We make 
the simplifying assumption that the decay amplitude fac- 
tors into currents C7 M (1 — 7s)b and S7^(l — 7s)c. Heavy 
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quark symmetry can then be used to express the 3-body 
double-charm decay amplitudes in terms of two functions 
Gi(g 2 ) and G 2 (q 2 ), where q 2 is the invariant mass of the 
hadrons produced by the S7 P (1 — 75)0 current [l^. For 
example, the amplitudes for decays into D°D*°K + and 
D* Q D Q K+ are 

A[B+ -> D Q D*°K + ] = -ide* ■ (V + v) 

-i(G 2 /m B K [«* ' k (V + v) u - «. ■ {V + v)k u 

-ie v ^(V + v)^ a k ], (12) 

A[B+ -> D*°D°K + } = i(G lVfl + G 2 kJm B )el 
x [(1 + v. ■ V)g^ - v»V v - ie^v^Vp] , (13) 

where k is the 4-momentum of the K + and V, v*, and v 
are the 4- velocities of the B + , D*° or D*°, and D° or D°, 
respectively. As a further simplification, we approximate 
G\ and G 2 by constants. The resulting expressions for 
the 3-body double-charm decay rates are 

T[B + -> D°D a K + ] = l(T 3 MeV 

x (178.9 |Gi| 2 + 51.8Re(G?G 2 ) + 4.37 |G 2 | 2 ) , (14) 

T\B+ -> D°D*°K + } = 10~ 3 MeV 

x (49.6 |Gi| 2 + 2.61 Re(G 1 G 2 ) + 3.49 |G 2 | 2 ) , (15) 

T[B+ -► D*°D°K + } = l(T 3 MeV 

x (52.5 |Gi| 2 + 1.87Re(GtG 2 ) + 2.31 |G 2 | 2 ) , (16) 

T[B+ -> D* a D*°K + ] = l(T 3 MeV 

x (221.5 |Gi| 2 + 74.8Rc(GtG 2 ) + 11.58 |G 2 | 2 ) (17) 

We obtain a good fit to the Babar branching fractions 
with Gi = 3.2 x 1CT 6 and G 2 = (-14.6 + 9.6i) x 10~ 6 . 
In the corner of phase space where the 4-velocities of D° 
and D*° are equal, the amplitudes (|12fl and (|13|l reduce 
to the form on the right side of JfjJ) with coefficients c\ = 
c 2 = —iGi/mB + iG 2 {nriB + rriD + mjj*)/m B . If X has 
charge conjugation C — +, the estimate (0 reduces to 



B[B+ -> XK+] w (2.6 x 10" 5 ) 



A? 



0.4 McV 



\ 1/2 

J .(18) 



Tx < 2300 keV. The sum of the widths for decay into 
D°D°ir° and D°D°j approaches r[L>*°] w 50 keV in the 

limit Ei, — > 0. The remaining partial widths scale as 

1/2 

EJ . Using a coupled-channel calculation in a model in 
which X mixes with J /if) p, the decay rate for J/ip 7r + 7r _ 
has been estimated to be 1290 keV for E b = 0.7 MeV [TT| . 
Thus it is at least plausible that J/ip 7r + 7r~ is one of the 
major decay modes. Other possible decay channels are 
T] c Trir, radiative transitions to charmonium states, and cc 
annihilation decays. 

We have calculated the decay rate for B + — * XK + 
and the differential decay rate for B + — * D a D*°K + near 
the D°D*° threshold under the assumption that X(3872) 
is a loosely-bound S-wave D D*°/D°D* molecule and 
that its production rate is dominated by the coalescense 
of charm mesons. Observation of a sharp peak in the 
D°D*° invariant mass distribution near threshold in the 
decay B + — > D°D*°K + would confirm the interpreta- 
tion of X as a D°D*° molecule. Our order-of-magnitude 
estimate of the branching fraction for B + — > XK + is 
compatible with observations if X(3872) has quantum 
numbers J PC = 1 ++ and if J/tp ir + ir~ is one of its ma- 
jor decay modes. 

This research was supported in part by the Department 
of Energy under grant DE-FG02-91-ER4069. 



If G = — , the branching fraction would be significantly 
smaller because of destructive interference between c\ 
and c 2 . We could get a more reliable result for the nu- 
merical prefactor in l|18[) by relaxing the factorization as- 
sumption and carrying out a Dalitz plot analysis of the 
3-body decays. Since the result depends quadratically on 
the ultraviolet cutoff Ai , the best we can do is obtain an 
order-of-magnitude estimate of the branching fraction by 
setting Ai « m^. 

The Belle collaboration measured the product of 
the branching fractions B[B + — * XK + ] and B[X — > 
J/tp 7T+TT-] to be (1.3 ± 0.3) x 10~ 5 1\. Our esti- 
mate of B[B + — > XK + ] is compatible with this result 
if J/tp 7r + 7r~ is one of the major decay modes of X. The [19] 
experimental upper bound on the width of A (3872) is 
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